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A LIAPUNOV  FUNCTION  AT,  FOR  A 
MATRIX  DIFFERENCE-DIFFERENTIAL  EQUATION 


Abstract:  A quadratic  positive  definite  functional  that  yields 

necessary  and  sufficient  conditions  for  the  asymptotic  stability 
of  the  solutions  of  the  matrix  difference-differential  erpiation 
x(t)  = Ax(t)  + Bx(t-x)  is  constructed  and  its  structure  analy/eil. 
This  functional,  a Liapunov  functional,  provides  the  best  possible 
estimate  for  the  rates  of  growth  or  decay  of  the  solutions  of 
this  equation.  The  functional  obtained,  and  its  method  of  con- 
struction, are  natural  generalizations  of  the  same  problem  for 
ordinary  differential  equations,  and  this  relationship  is 
emphasized.  An  example  illustrates  the  applicability  of  tfie 
results  obtained. 


1 . Introduction. 

In  this  paper  we  construct  a functional  for  the  determination 
of  the  asymptotic  behavior  of  the  solutions  of  the  linear  autonomous 
matrix  difference-differential  equation  with  one  delay 


x(t)  = Ax(t)  + Bx(t-T),  t > 0, 


where  x(t)  is  an  n-vector  function  of  time,  A and  B are  con-^i » 
stant  n X n matrices  and  x > 0. 

In  a recent  paper  [7]  such  a functional  was  constructed 
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for  the  special  case  when  x(t)  is  a scalar  function;  there,  the 
Liapunov  functional  was  obtained  as  the  limit,  in  an  appropriate  “•f’ 


n □ 


2 


► 
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sense,  of  a Liapunov  function  constructed  by  well-known  methods 
for  a difference  equation  approximation  of  the  original  functional 
equation.  In  this  paper,  we  extend  the  results  in  [7]  so  that 
they  are  applicable  to  the  general  matrix  case;  we  shov/  that  the 
general  form  of  the  functional  in  the  matrix  case  is  a natural 
generaliiiation  of  the  one  obtained  in  [7]  and  tliat,  with  appropriate 
modifications,  all  of  the  arguments  presented  in  that  paper  are 
applicable  here.  Essential  to  this  generalization  was  the  study, 
presented  in  [2],  of  the  existence,  uniqueness  and  structure  of 
the  solutions  of  a special  functional  differential  equation  which 
is  intimately  connected  with  the  construction  of  the  Liapunov 
functional  presented  here. 

The  Liapunov  functionals  obtained,  and  their  method  of  con- 
struction, are  very  natural  extensions  of  the  methods  used  for 
the  construction  of  Liapunov  functions  for  linear  systems  of 
ordinary  differential  equations.  We  attempt  to  bring  this  relation- 
ship into  evidence. 

A simple  example  illustrates  the  applicability  of  the  results 
obtained,  by  giving  a detailed  construction  of  a functional  for  a 
2 X 2 system. 

The  problem  of  the  construction  of  Liapunov  functionals  for 
equations  of  this  type  has  been  previously  considered  by  Repin  [8] , 
Datko  [3]  and  Hale  [5,6].  This  study  is  in  the  spirit  of  the 
previous  ones,  but  is  more  specific.  The  functionals  obtained  here 
are  more  general  than  the  ones  used  by  Hale,  and  they  yield  not 
only  sufficient  but  necessary  conditions  as  well  for  asymptotic 
stability.  The  paper  of  Repin  does  not  seem  to  us  to  be  correct. 


J 
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through  a simple  mistake  with  profound  consequences.  Patko  [3J 
has  not  pursued  the  problem  of  construction  to  the  extent  presented 
here,  but  his  ideas  have  had  a deep  effect  on  our  approach. 


2 . The  Difference-Differential  Equation. 

Denote  by  ( [a  ,b]  , the  space  of  Lebesgue  square 

integrable  functions  defined  on  [a,b]  with  values  in  , and  for 

a fixed  t > 0 consider  the  Hilbert  space  ^ x L~  ( [-t  , 0]  , 

-0 


with  inner  product  <u^^,U2>  = Vj^V2 


-T 


4)^  (9 ) (})2  { 0 ) d0  , where 


u.  = (v.,(b.)  € and  the  naturally  induced  norm  ||(v,4i)JI 
111 


V V + 


(j)'^(G)  4)  (0)dO.  With  x:  [-x,''’)  >■  , for  t ;;;  0 wo  denote 


-X 


by  X the  function  x : [~x,0]  ->■  , where  x (0)  = x(t+9). 

U t t* 


Consider  the  matrix  difference-differential  equation 


x(t)  = Ax(t)  + Bx(t-x),  t > 0,  (2.1) 

where  A,B  are  n x n matrices,  x(t)  is  an  n-vector  and  x > 0, 
together  with  the  initial  conditions 


Xp(0)  = C,  Xp  = 4).  (2.2) 

where  (C»4>)  € dl. 

A solution  of  this  initial  value  problem  is  a function 
X e L2 ( [-T , t] , for  each  t > 0,  such  that  x is  absolutely 
continuous  for  t > 0,  it  satisfies  (2.1)  a.e.  on  [0,t]  and 
x(0)  = x(0)  ==  4>(0)  a.e.  for  0 G[-t,0].  It  is  known  [1,6] 

that  (2.1)-(2.2)  has  a unique  solution,  defined  on  [-x,»),  which 
depends  continuously  on  the  initial  data  in  the  norm  of  tH. 
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The  initial  value  problem  (2.1) -(2. 2)  can  be  rewritten  as 


dt 


'x^(O)' 

= V 

x^(0)' 

X. 

[ t J 

t 

(2.3) 


(Xq(0)  ,Xq)  = (C,<P)  e (2.4) 

where 


x^(0) 

Ax^(O)  + Bx^(-x) 

II 

9x  (0) 

l"t  , 

1-T15— i « ; »J 

and  this  operator  has  a domain,  dense  in  defined  by 


^(J^)  = {(^,(p)  e <P  is  A.C.  in  [-1,0], 

({)'  G L2[-t,0]  , (^(0)  = r,}. 

The  operator  ^ is  the  generator  of  the  CQ-semigroup  .9^'(t)  , where 
y\t)  : ii:  ^ is  given  by  ^(t)(C,(p)  = (x(t),x^),  (x(t),x^)  the 

solution  pair  of  (2.3),  (2.4). 

It  is  known  [1,6]  that  there  is  a constant  y such  that  the 
spectrum  of  lies  xn  the  left-half  plane  Re(X)  < Y»  and  that 

for  every  G > 0 there  exist  a constant  K > 1 such  that 

I I ^ (2.6) 

The  spectrum  of  ja^  consists  of  those  complex  X that  satisfy  the 
characteristic  equation 
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det[AI-A-Be  = 0. 


(2.7) 


Finally, [5,6],  a useful  representation  of  the  solutions  of  (2.1) 
is  given  for  every  t,  T > 0,  by 


x^_^^(0)  = S(T)x^(0)  + 


-T 


S (T-a-r)  Bx^  (oi)da. 


(2.8) 


where  the  matrix  S is  the  solution  of  the  matrix  initial  value 
problem 


= S (t) A + S (t-T)B, 


S(0)  = I,  S(t)  = 0 for  t < 0. 


(2.9) 


We  now  turn  to  the  construction  of  a Liapunov  functional. 

3 . A Quadratic  Functional. 

Associated  with  the  functional  differential  equation  (2.1), 
or  (2.3)-(2.4),  and  motivated  by  the  results  in  [7],  we  wish  to 
consider  the  real  symmetric  quadratic  form  on  where 

([.A)  e 


V(C»<]>)  = + e*^^  f (|)'^(0)Re^'^°<))  (O)d0  + 

J -T 

+ f/Q(0)^  + 25^  [ Q(a+T)e‘^^“'^^^B(J)(a)da  (3.1) 

* -T 

+ 2 [ [ (})'^(a)B'^Q(6-a)e‘^  ^“■^^'^^■^^B())(e)d6da, 

''  -T''  a 

where  6 is  a real  number,  M,R  are  constant  n x n real  positive 


A 


-I 
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definite  matrices  and  Q(a)  is  a continuously  differentiable  matrix 
that  is  assumed  to  satisfy  the  initial  value  problem  for  the 
functional  differential  equation 

Q'  (a)  = (a'^+6I)Q(«)  + e'^’^sV  (x-a)  , 0 < a < x (3.2) 

Q(0)  = Q(O)'^  = Qq,  (3.3) 

where  the  superscript  T denotes  transpose  and  is  a symmetric 

but  otherwise  arbitrary  matrix. 

Evaluation  of  the  Liapunov  functional  (3.1)  aloinj  tlic  solu- 
tions of  (2. 3) -(2. 4)  yields  a function  of  time,  which  we  denote  by 
V(t)  = V(x^(0),x^).  This  function  of  time  is  differentiable  along 
the  solutions,  and  a laborious  but  straightforward  computation,  which 
mal<es  use  of  (3.2)  and  (3.3),  shows  that  the  derivative  of  this 
function  along  the  solutions  is  given  by 

V(t)  = dt  '^’t^  " -26V(x^(0)  ,x^)  + (3.4) 

+ x^(0)  [ (A^+6I)Q(0)  + Q(0)  (A+6I)  + (A^+6I)M  + M(A+6D  + 

+ e^^B'^Q'^(T)  + e'^^Q(x)B  + 2e‘^'^R]  x^  (0) 

- e'^^h  (0)  ,x^  (t)  ) , 

where 


R MB 

X.  (0) 

h(x^.(0)  ,x^(-t))  = [x^(0)  ,-e"'^'^xJ{-T)] 

T 

B M R 

c 

-6r  , . 

[-e  x^.(-T^ 

(3.5) 

t 

7 

It  is  our  purpose  to  show  that,  through  the  functional 
(3.1)  and  its  derivative  (3.4)  it  is  possible  to  estimate  the 
asymptotic  behavior  of  the  rate  or  growth  on  decay  of  the  solutions 
for  our  original  functional  differential  equation  (2. 3) -(2. 4). 

For  this  purpose,  appropriate  choices  must  be  made  for  the  positive 
definite  matrices  M and  R and  of  the  constant  6,  and  of  the 
matrix  Q(a)  determined  by  (3.2)  and  (3.3). 

To  be  more  specific,  denoting  by  y = max(Re  X|  det[XI-A- 
Be  ] = O),  we  wish  to  show  tliat  for  every  (E  > 0 and 
-6  = Y + 2G  it  is  possible  to  choose  matrices  M,R  and  Q(a) 
satisfying  (3.2)  and  (3.3)  so  that  there  exist  positive  constants 
c^,C2  such  that 

< V(5,4>)  < c I I (f;,(|.)  I 1^^  , (3.6) 

and 


V(C,<}>)  < -26V(5,4>)  . 


(3.7) 


1/2 

These  last  relationships  imply  that  {V(^,<f>)}  = | | u ( ■ , i,*' ) | | is  a 

norm  equivalent  to  the  original  norm  on  and  that,  in  this  norm 


(x^(0) ,x^)  I 1^  < I |Xq(0) ,Xq|  I 


e-5^ 


(3.8) 


whereas,  in  the  original  norm  on  , 


rc--vi/2 


< 


(3.9) 
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These  estimates  are  precisely  those  stated  in  (2.8)  and  are  the 

best  possible  ones.  It  should  bo  noted  tliat  the  lujrm  itulucfnl  l)y 

the  square  root  of  the  Liapunov  functional  is  the  best  possible 

one  in  the  sense  that  it  yields  (2.6)  with  K = 1.  Moreover,  if 

Y <0,  the  above  remarks  show  that  the  Liapunov  functional  (3.1) 

yields  uniform  exponential  asymptotic  stability. 

Consider,  first  of  all,  the  appropriate  choice  for  the 

matrix  Q(n),  a solution  of  (3.2),  (3.3)  for  0 v u . 1 n [ .’ | 

it  was  shown  that  equation  (3.2)  with  initial  conditions  (3.3) 

has  a unique  solution;  moreover,  it  was  shown  that  the  linear 

2 

vector  space  of  all  solutions  of  (3.2)  has  dimension  n . Indeed, 
introducing  the  notation 


Q (a)  = (q^  j (cx)  ) 


^1* (a) 


[q^j^  (a)  , . . . ,q*j^  («)  1 , 


where  q.*(a)  and  q* . (a)  are,  respectively,  the 

1 * j 

t hi  2 

the  j column  of  Q(a),  defining  the  n -vector 


. th 
1 


row  and 


T 

q(a)  = (a)  , . . . (a)  ) , 

T 

and  defining  the  matrix  R(a)  = Q (i-a),  equation  (3.2)  reduces, 

through  the  use  of  the  Kronocker  (or  direct)  product  of  two 

2 

matrices,  to  the  2n  system  of  ordinary  differential  equations 


q(a) 

(A+61)'^  0 I 

0 I 

q (a) 

da 

r(a)_ 

-1  0 

-I  0 (A+ai)"^ 

r (a)_ 

with  the  auxiliary  condition 


(3.10) 


f'inj'''  = 


where  F is  an  arbitrary  n x n matrix.  The  solution  of  (3.10), 

2 

subject  to  (3.11) , yields  solutions  Q(a)  with  n parameters 

that  can  always  be  chosen  to  satisfy  any  arbitrary  condition  of  the 

T 

typo  Q(0)  = Q(0)  = Qq.  In  [2]  the  structure  of  these  solutions 

and  simple  methods  of  computation  that  take  advantage  of  the 
particular  structure  of  the  equation  are  presented. 

Associated  with  equation  (3. 2) -(3. 3)  is  an  integral,  v;hoso 
structure  is  motivated  by  a similar  integral  for  ordinary  differ- 
ential equation.  Indeed,  let  W be  a symmetric  matrix  and  let 
S(t)  be  the  solution  of  equation  (2.9).  Consider  the  expression 


Q(a)  = (u)  e'^^^S  (u-a)  e"^  '^^du. 

0 


(3. 12) 


Since,  for  every  €:  > 0 , | |S(t)  | | < for  some  K > 1 and 

since  5 = -y  - 2G,  it  follows  that  this  integral  converges.  More- 
over, it  immediately  follows  from  the  definition  of  Q(a)  that 
Q(a)  = o'^(-a),  5(0)  = q'^(O)  and,  since  S(t)  satisfies  (2.9), 
that 


5'  (a)  = -Q(n)(A+6l)  - 5(a+T)e'^'^B  - (a ) e'^'^'w. 


Since  Q'  (a)  = ^ [o'^(-a)]  for  a 0,  we  also  obtain  that 


5'  (a)  = (A+6I)'^Q(a)  + (r-a)  + s'^(-a)e^"w. 


From  these  relationships  it  follows,  since, for  a > 0,  S(-a)  = 0, 
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that  Q(a)  satisfies 


Q'(a)  = (A'^+6I)Q(a)  + (i-a)  , 0 < a - r 

#CO 

q'^(O)  = Q(0)  = s'^(u)o'^%S(u)e'^^du 

J 0 

and,  moreover,  given  the  continuity  of  Q(a) , that 

Q'(0)  + Q''^(0)  = (a'^+6I)Q(0)  + Q(0)  (A+6I)  + b'^g'^'^q'^  ( ; ) 
= -s'^(0)W  = -W. 


(3.13) 


+ QCOe'^'^B 
(3.14) 


These  observations , and  the  uniqueness  of  the  solutions  of 
(3 . 2)  - ( 3 . 3) , show  that  Q(n)  as  defined  by  (3.12)  is  I he  uniijiu' 
solution  of  (3.2)-(3.3)  with  the  initial  conditions  proscribed  by 
the  second  equation  of  (3.13). 

It  is  transparent  that  to  each  constant  symmetric  positive 
definite  matrix  W there  corresponds  a Q(n)  given  by  (3.12) 
and  that  this  is  the  unique  solution  of  (3. 2) -(3. 3)  with  initial 


condition  Q(0)  = Q(0)  = 


(u)  (u)  e'^'‘*du  which  is  positive 


definite.  Conversely,  to  each  Q(0)  = Q (0)  = Q^,  the  unique 
solution  of  (3.2)-(3.3)  yields  a differentiable  matrix  function 
Q(a)  that,  through  (3.14)  defines  a unique  symmetric  matrix  W 
so  that,  for  this  W,  (3.12)  yields  the  unique  solution  to  (3.2)- 
(3.3)  with  the  prescribed  initial  conditions.  These  remarlcs  show 
that  the  map  W ->  Q(0),  defined  by 


(u)  e*^  ^WS  (u)  e'^^du, 


0 


Q(0) 


0( 


as  a map  on  the  space  of  n x n symmetric  matrices  is  one-to-one, 
onto,  and  it  maps  positive  definite  matrices  W into  positive 
definite  matrices  Q(0). 

With  this  particular  characterization  of  the  matrix  function 
(a)  it  is  possible  to  bring  into  evidence  the  particular  structure 
of  the  Liapunov  functional  (3.1).  Indeed,  substitution  ul  (3.1.’.)  into 
(3.1)  for  Q(a)  yields,  after  some  rearrangements  and  interchanges 
of  integrals,  that 


T <5t/T  260 

V(^,(p)  = C + e ^ rf,M0)Re  <{i(0)d0 


-T 


{S(u)f  + 


0 

0 

J 


S (u-ck-t ) B(J>  (a ) d K ^ B ( 11 ) ' 


(3.1 1 


S (u-ot-T  ) B<))  (a)  da  }du , 


or,  evaluating  this  functional  on  tl\e  solutions  of  our  equation, 
and  using  (2.8)  , 


V(x^(0) 


,Xt)  = 


xj’(0)tlx^(0)  + e'^'^ 


-T 


T 2 (S  0 

x^(0)Re  x^.  (0)d0 


(3.1G) 


Similarly,  in  this  notation,  equation  (3.4)  becomes 


V(x^(0) ,x^)  = 


-26V(x^(0)  ,x^)  +x^(0)[-W+  (a'^+6I)M  + M(A+6I)  + 


+ 2e'^^R]x^(0)  + [x^(0),-e  '■’‘xJ.(-t)] 


-6t  T, 


R MB 


T 

B M R 


x^(0) 

-6t  T,  , 
-e  x^(-T) 


(3.17) 
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Given  the  nonnegative  nature  of  the  last  term  in  (3.16)  it 


follows  that 


min  (M)  ,G  ^ ^ ^ ^niin  ^ ^ ^ ^ ^ ^ 


yielding  a value  of  c^,  for  equation  (3.6),  given  bv 


c = min(X  (M),e  . (R)). 

1 min  min 


(3.18) 


Similarly,  since,  from  (2.6), 


x.,(0)||  _ < II  (x„  (0)  ,x,J  I (_  < I (Xq(0)  ,Xj^) 


, n - ' ' u u ' ~ 


it  follows  that 


yielding,  for  equation  (3.6),  a value  of  given  by 


I r I ^ ( W ) « 

r ^ /.«x  <5  T , max  „2 


(3.19) 


It  remains  to  be  shown  that  equation  (3.7)  holds.  For  this 
purpose  inspection  of  equation  (3.17)  indicates  that  it  is  necessary 
and  sufficient  to  show  that  it  is  possible,  by  appropriate  clioicos 
of  positive  definite  matrices  VJ,R  and  M,  to  have  the  matrix 
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G = 


W - (A+5I)'^M  - M(A+6I)  - Re'^’' 


(5t 

B Me 


MBe 


6t 


Re 


6 T 


(3.20) 


positive  scmidof inite , But  this  is  always  possible;  ,i 

particularly  simple  choice  of  such  matrices  is  to  lot  W,R  and 
M to  be  nonnegative  multiples  of  the  identity  matrix,  i.o. 

M = I,  R = and  W = With  >>  k^^  >>  1,  the  matrix  G 

is  clearly  positive  definite,  hence  its  determinant  is  positiv^e, 
and  the  matrix  will  become  semidefinite  only  when  the  determinant 
vanishes.  But,  for  the  above  choice  of  matrices  wc  have 


det{G)  = det(kj^Ie’^'^)dot  [k^I  - (A+A*^)  - 261  - e'^^  (k^^I  + ^ b'^B)  J ; 

R 


Letting  k = (B^B)  and 


R max 


k^^  = max{0,  (A+a"^)  + 26  + 2o'^’^  } , 


max 


max 


then  (3.20)  is  positive  semidefinite  and  tiiorefore,  for  those 
choices  of  M,R  and  W, 


V(x^(0),x^)  < -26V(x^.  (0)  ,x^)  . 


(3.21) 


It  is  to  be  remarked  that,  with  these  choices,  the  Liapunov  func- 
tional constructed  is  particularly  simple  and  reduces  to 
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V(x^(0),x^)  = 
+ 


This  functional  is  a direct  generalization  of  that  used  in  [7],  and 
generalizes  those  used  by  Hale  [5,6]  and  Datko  [3].  Wo  recapitulate 
the  above  results  in  the  form  of  a 

Theorem;  Given  the  retarded  equation  x(t)  = Ax{t)  + Bx(t-r) 
and  the  Liapunov  functional  V given  by  equation  (3.1),  if 
Y = max {Re  X|  det [ A I-A-Be"^^ ] = 0}  and  G > 0,  then  there 
exist  constant  positive  definite  matrices  M and  R and 
a differentiable  matrix  Q(n),  0 < a < t with  Q(0)  = Q(0)"^ 
such  that  the  functional  V is  positive  definite,  bounded 
above,  and  V < 2(y+G)V. 

Of  course,  if  Y 0,  the  above  result  implies  exponential 
asymptotic  stability;  moreover,  the  rate  of  decay  is  precisely  the 
expected  one. 

4.  A Comparison  with  Ordinary  Differential  Equations. 

It  is  our  object,  in  this  section,  to  point  out  the  intimate 
relationship  between  the  results  obtained  for  our  functional  equation 
and  the  classical  results  for  the  construction  of  Liapunov  functions 
for  ordinary  differential  equations. 


x^(0)x.  (0)  + e'^'^  A (b'^B) 
t t max 


x|^(0)x^(0)e^'^®d0  + 


max{0,A  (A+A*^)  + 26  + 2e"'^/A  (b'^’b)}.  (3.22) 

' max  max 


' T 


x^  (0)x  (O)e^'^^ds. 

t+s  t+s 
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Recall,  14],  that  given  the  system  of  ordinary  differential 
equations 


x(t)  = Cx  (t)  , (4.1) 

where  C is  an  n x n matrix,  a Liapunov  function  for  this 
system  can  be  always  taken  as  the  quadratic  form 

V(x)  = x'^Px,  (4.2) 

where  f is  a positive  definite  matrix.  Moreover,  if 
Y = max{Re  A|  det(AI-C]  = 0}  and  if,  for  G > 0,  6 = y + 2G, 
then  given  an  arbitrary  positive  definite  matrix  W,  the  algebraic 
equation 


(0+61)"^?  + P(C+6I)  = -W, 


(4.3) 


has  a unique  solution  P which  is  positive  definite.  This  matrix  p, 
if  used  in  (4.2)  along  the  solutions  of  the  differential  equation 
(4.1)  yields,  upon  differentiation,  that 


^(x(t))  = -26V(x(t))  - x'^Wx  < -26V{x(t)). 


(4.4) 


Furthermore,  the  unique  positive  definite  solution  P of  (4.3)  can 
be  obtained  as  the  integral 


= f: 


P 


(4.5) 
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Let  us  now  bring  into  evidence  the  relntionsl\ip  between  the 
results  obtained  in  the  previous  section  and  these  results.  The 
functional  dif  f erent  j.al  equation  under  eonr;  i derat  i on  in 


x(t)  = Ax(t)  + Bx(t-x),  t > 0 . 


(4.1'  ) 


“At 

Once  again,  we  assume  tliat  y = max{Re  X|  det[M-A-Be  } = 0} 
and,  for  G > 0,  let  -6  = y + 2€.  The  Liapunov  functionn]  is 
then  of  the  form 


V(x^(0),x^)  = xJ|(0)Mx^{0)  + e'^'^  xj  ( 0 ) Re^ '^°x^  { 0 ) dO  (4.2') 

-T 


+ x^(0)Q(0)x^(0)  + 2x'^(0)  Q(a+t)e'^ 


Bx^ ( T ) do  + 


0 rO 


x^(a)B'^Q(3-«)e'^  ^“■^^■^^■^^Bx^(p)dpda. 


The  choice  of  the  positive  definite  symmetric  matrices  M and  R 

in  this  expression  is  rather  arbitrary.  Their  purpose  is  to  insure 

the  strict  positive  definiteness  of  the  functional  on  the  Hilbert 

space  ; it  should  be  noted  that  if  M = R = 0,  the  functional 

(4.2')  is  positive,  but  does  not  satisfy  a relationship  of  the  type 

V(x.  (0),x.)  > c,  I I (x.  (0) ,x. ) I I for  c,  > 0.  The  requirement  that 
u t “ X tl  t ^ X 

the  matrix  G of  equation  (3.20)  be  positive  seraidefinite  is  always 
satisfied  for  M = R = 0 and  W positive  semidef inite.  Given  an 
arbitrary  positive  definite  matrix  W,  it  is  always  possible  to 
select  positive  definite  matrices  M and  R so  that  G is  positive 
semidef inite. 
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The  choice  of  the  continuously  differentiable  matrix  0(,y), 
0 < (X  < r is  critical.  It  must  satisfy,  given  an  arbitrary 
positive  definite  matrix  W,  the  functional  equation 


Q'(a)  = (A'^+6l)Q(a)  + (t-u)  , 0 < a < i, 

Q(0)  = q'^(O)  , 


(4.3') 


with  the  condition 


(a’^+6I)Q(0)  + Q(0)  (A+6I)  + + Q(x)e‘^’^B  = -W.  (4.3") 


In  the  previous  section  it  was  shown  that  such  a Q(a)  always 
exists  and  is  unique. 

With  such  a choice  of  Q(a) , then  one  obtains  that,  along 
the  solutions  of  the  functional  differential  equation  (4,1') 

V(x^(0),x^)  <-26V(x^(0) ,x^) . (4.4') 

Moreover,  such  a matrix  Q(a)  exists,  is  unique,  and  a representation 
of  it  is  given  by  the  integral 

Q(n)  = f s'^(u)e'^“ws  (u-n)e‘^  ^^-''‘^du,  (4.3') 

J 0 

c 

where  S(t)  is  the  solution  of  equation  (2.9). 

The  strong  relationship  between  the  unprimed  and  primed 
equations  is  now  clear.  Indeed,  note  that  for  t = 0,  the  matrix 
C in  (4.1)  becomes  A + B and  the  matrix  V in  (4.3)  becomes 
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Q(0)  + M and  all  the  primed  equations  become  the  unprimod  ones. 

Equations  (4.3’)  and  (4.3")  are  of  a mucli  more  complex 

nature  than  the  familiar  algebraic  equation  (4.3).  However,  in 

spite  of  its  appearance,  the  linear  vector  space  of  the  solutions 

of  (4.3')  is  not  infinite  dimensional  but,  as  was  pointed  out  in 

2 

the  previous  section,  has  dimension  n . Hence,  although  the 
problem  of  construction  of  the  Liapunov  functional  for  the  func- 
tional differential  equation  does  not  reduce  to  the  solution  of 
an  algebraic  equation  such  as  (4.3);  it  reduces  fo  the  soiulion  ot 

a linear  differential  equation,  equation  (3.10),  which  is 
accomplished  without  difficulty. 

5 . An  Example. 

To  illustrate  the  application  of  the  method  developed  hero, 
we  wish  to  briefly  outline  the  construction  of  an  appropriate 
Liapunov  functional  for  a simple  two  dimensional  example.  For 
this  purpose,  consider  the  system 


-5  0~ 

~4  2 

x(t)  = Ax(t)  + Bx(t-r)  = 

x(t)  + 

_ 0 2_ 

_0  0_ 

where  t = 6.  In  the  notation  of  Section  3,  the  value  ~ 
is  selected  as  is  indicated  by  a simple  examination  of  the 
characteristic  equation  of  the  system.  For  the  construction  of 
the  Liapunov  functional  (3.1),  appropriate  positive  definite 
matrices  M and  R must  be  selected  and  a matrix  function  Q(a) 
that  satisfies  (3.2),  (3.3)  and  (3.14)  must  be  determined  given  a 
preselected  positive  definite  matrix  W. 
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As  indicated  in  Section  3,  an  appropriate  choice  for  the 
matrices  M,R  and  W is  ctivcn  by 


M = I, 


(5.2: 


W = max{0,A  (A+A^)  + 26  + 2e‘^\A  (B^n))l  = 7,24r,I  = k,,I. 

' max  max  w 


All  that  remains  is  to  determine  the  unique  matrix  function 
Q(a)  that  satisfies  (3.2),  (3.3)  and  (3.14).  The  general  form 

of  this  matrix  function  is  determined  by  obtaining  the  solutions 
of  the  differential  equation  (3.10)  following  the  methods  indicated 
in  [21 ; in  a straightforward  although  somewhat  lengthy  computation 
one  obtains  that  the  general  form  of  this  matrix 


Q (a)  = 


q3^l(a) 

q2i(a) 

^22 

, 0 < a < T , 


(5.3) 


is  of  the  form 


qj^^(a)  = B, 


qi2(ot) 


^/fW)  (.a-T/2)  + (“-T/2) 


+ 34  e 


(2-6) («-t/2) 
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(-5+6)  (u-i/2) 

I'"  ' '^2 


q2^(a)  = +3 


H22^  (5 ) e*' 


£T'6j  (.t-i/2)  , 


-/fT6) (a-T/2) 


+ 3^g2('5)c 


-(2-6) («-t/2) 


and 


*^22^”^  ^2 


-/f  ( 6 ) (a-x/2)  , ,5.. 

nT-,(6)e  + (6)  H22  J ^ 


-/fT6)  (cx 


'22 


+ («-T/2)  ^ 3^nii(6) 


(2-6) (a-x/2)  ^ 


+ (6)  e 


-(2-6) (a-x/2) 


In  these  expressions  the  indicated  functions  of  6 are, 
given  by 


f (6) 


(5-6)^  - 16e^'^\ 


g]^('5 ) 


_ /TW)  + (5-6) 


4e 


6x 


92(^)  = 


(6) 


ri2i(6) 


7-26 

7- 26 

8- 46 


2e 


6x 


[2-6-/FTT)  ]gj^(6) 

2e'^'^gj^(6) 


2-6+/f (6) 


(5.  -1) 


-t/2)  ^ 


respectively , 


(5.5) 


i 
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The  constants  ^4  determined  from  the 

four  equations  obtained  from  the  relation  Q(0)  = Q(0)^  (i.e., 

qi2(0)  = from  equation  (3.14)  (i.e,,  O'(0)  4-  Q'(O)’^ 

= -w  = -V). 

Solution  of  these  four  linear  equations  for  the  indicated 
value  of  6 = ^ yields  the  values 

= -O.llS.lO'^kj^,  ^2  = 0. 321. 10~^k^^, 

(5.6) 

^3  = 0.242.  lO'^k^^,  = -0. 622. 10"^k^^. 

With  these  values,  it  is  easily  computed  that  Q(0)  = Q(O)'^,  and 
that  this  matrix  is  positive  definite. 
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